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LOS  ANGELES  AND  LONG  BEACH  HARBORS.  CALIFORNIA 
MODEL  ENHANCEMENT  PROGRAM 

NUMERICAL  WATER  QUALITY  MODEL  STUDY  OF  HARBOR  ENHANCEMENTS 


Introduction 


Background 

1.  Los  Angeles  and  Long  Beach  Harbors  (San  Pedro  Bay),  California,  are 
located  adjacent  to  each  other  and  share  a  common  breakwater  system  (Fig¬ 
ure  1).  The  Ports  of  Los  Angeles  and  Long  Beach  have  undertaken  a  long-range 
cooperative  planning  effort  known  as  the  2020  Plan.  Incorporated  in  the  plan 
are  harbor  enhancements  that  include  dredging  and  new  landfill. 

Problem  identification 

2.  Harbor  enhancements  may  affect  water  quality  in  the  study  area  by 
changing  the  tidal  circulation  and  flushing  patterns.  The  major  water  quality 
concern  was  the  dissolved  oxygen  (DO)  resource.  Channel  deepening  introduces 
the  possibility  of  DO  stratification  where  well-oxygenated  surface  waters 
overlie  oxygen-depressed  bottom  waters. 

Objective 

3.  The  purpose  of  this  study  was  to  compare  the  flushing  and  DO 
•esources  of  existing  and  plan  conditions  through  numerical  model  simulations. 

A  three-dimensional  (3-D)  water  quality  model  (WQM)  was  used  to  address  possi¬ 
ble  vertical  stratification.  Hydrodynamic  output  from  the  3-D  hydrodynamic 
model  (HM)  CH3D  was  used  as  input  by  the  WQM  to  transport  water  quality 
constituents . 

Report  organization 

4.  The  succeeding  sections  of  the  report  describe  modifications  made  to 
the  WQM  code  and  the  procedure  used  to  interface  the  WQM  and  CH3D:  the  exer¬ 
cises  used  to  ensure  that  the  transport  properties  of  CH3D  were  maintained 
with  the  interfacing;  and  the  flusliing  and  water  quality  studies.  The  water 
quality  kinetic  routines  are  detailed  in  Appendix  A. 
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Figure  1.  Location  of  channels  and  basins  in  Los  Angeles  and 
Long  Beach  Harbors,  California 


Water  Quality  Model 

5.  The  water  quality  model  selected  for  this  study  was  a  modification 
of  the  US  Environmental  Protection  Agency  WASP  code  (Ambrose,  Vandergrift,  and 
Wool  1986).  Major  adaptations  included  (a)  improved  advective  and  diffusive 
transport  schemes,  (b)  provisions  for  the  input  and  processing  of  CH3D  3-D 
hydrodynamic  data,  and  (c)  implementation  of  kinetic  routines  specific  to  the 
San  Pedro  Bay  application. 

Advective  and  diffu- 
sive  transport  schemes 

6.  The  original  WASP  model  formulation  considers  both  advective  and 
diffusive  transport  through  faces  of  adjacent  cells  of  arbitrary  size,  shape, 
and  distribution.  For  the  present  study,  horizontal  flows  were  distinguished 
from  vertical  flows,  and  improved  transport  schemes  were  implemented. 

7.  The  advection  scheme  used  in  the  original  WASP  provided  three  advec- 
tion  options:  upwind  difference,  central  difference,  or  a  linearly  weighted 
mixed  upwind/central  difference  scheme.  These  options  were  unsatisfactory 
because  they  either  introduced  numerical  diffusion  to  the  extent  that  physical 
diffusive  and  dispersive  transport  were  entirely  masked  or  required  exces¬ 
sively  large  values  of  physical  diffusion  to  ensure  numerical  stability.  To 
improve  the  horizontal  advective  transport  scheme  used  in  WASP,  a  modified 
version  of  the  QUICKEST  (Quadratic  Upstream  Interpolation  for  Convective  Kine¬ 
matics  with  Estimated  Streaming  Terms)  scheme  (Leonard  1979,  Hall  and  Chapman 
1985)  was  implemented.* 

8.  The  vertical -  stretching  relationship  in  the  hydrodynamic  model 
accounts  for  much  of  the  vertical  advection  through  layer  thickness  expansion 
or  contraction.  As  a  consequence,  the  vertical  velocities  were  small.  Small 
vertical  velocities  result  in  small  numerical  diffusion  with  the  upwind 
scheme,  and  large  physical  diffusion  is  not  required  to  maintain  stability 
with  the  central  difference  scheme.  However,  the  time  step  was  limited  due  to 
the  explicit  vertical  diffusive  transport  schemes  used  in  the  original  WASP. 
Therefore,  an  implicit  vertical  advective  and  diffusive  transport  scheme  was 


*  Ray  Chapman  and  Associates.  1988.  "Analysis  and  Improvement  of  the  Numer¬ 
ical  and  Physical  Mixing  Characteristics  of  the  WASP  Box  Model,"  report  pre¬ 
pared  for  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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implemented.  Central  differences  were  used  for  both  the  advection  and  diffu¬ 
sion  terms. 

Model  linkage 

9.  The  hydrodynamic  and  water  quality  models  were  linked  by  spatially 
and  temporally  averaging  CH3D  output  to  drive  the  WQM.  The  hydrodynamic  model 
used  extensive  spatial  resolution  to  resolve  important  scales  and  to  minimize 
the  need  for  parameterization.  The  CH3D  spatial  resolution  was  of  the  order 
of  100  m  and  required  a  time  step  of  60  sec  for  stability.  In  contrast,  the 
WQM  has  characteristic  time  scales  on  the  order  of  hours,  which  are  determined 
by  the  kinetic  rate  coefficients.  The  desired  analyses  of  water  quality  allow 
a  spatial  resolution  an  order  of  magnitude  less  than  that  used  by  CH3D. 
Reductions  in  spatial  and  temporal  resolution  greatly  reduced  computer  time 
required  for  the  WQM. 

10.  Plate  1  is  a  schematic  of  the  CH3D  grid  used  for  the  existing  con¬ 
ditions,  The  CH3D  grid  consisted  of  three  vertical  layers.  Existing  condi¬ 
tions  are  referenced  in  this  report  as  "Base."  Since  the  WQ.M  requires  less 
spatial  resolution  than  CH3D,  spatial  characteristics  of  the  hydrodynamic  grid 
such  as  volume,  cell  surface  areas,  cell  facial  areas,  and  cell  lengths  were 
summed,  resulting  in  an  "overlaid"  WQM  grid.  Plate  2  represents  the  resulting 
WQM  grid  overlying  the  CH3D  grid  for  the  Base  condition.  The  WQM  grid  main¬ 
tained  the  same  vertical  resolution  as  the  CH3D  grid  of  three  vertical  layers. 
Plate  3  represents  the  overlaid  WQM  grid  for  Scheme  B,  Phase  I,  of  the 

2020  Plan  enhancement,  which  will  be  referred  to  as  "Scheme  B." 

11.  Coding  was  appended  onto  CH3D  to  initially  write  to  an  output  file 
the  time  -  invar iant  data  that  included  cell  surface  areas  and  horizontal  cell 
lengths  and  to  temporally  average  horizontal  flows  and  facial  areas  and  write 
this  information  and  volumes  to  the  output  file.  The  hydrodynamic  data  were 
averaged  over  1-hr  intervals.  Averaging  over  shorter  intervals  did  not  alter 
WQM  transport  results.  A  mapping  file  was  prepared  that  was  used  by  the  CH3D 
appended  code  to  relate  the  hydrodynamic  and  water  quality  cells. 

12.  Table  1  displays  the  number  of  computational  cells  for  both  CH3D 
and  the  WQM.  Overlaying  the  WQM  grid  decreased  the  number  of  computational 
cells  by  94  percent.  The  larger  WQM  cell  sizes  allowed  an  increase  of  the 
time  step  used  from  60  to  900  sec  while  maintaining  computational  stability. 
Reductions  in  spatial  resolution  greatly  reduced  the  computer  time  required 
for  the  WQM. 
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Table  1 

Comparison  of  the  Number  of  Computational  Cells  in  CH3D  and  the  WQM 


CH3D 

WQM 

Plan 

Layers 

Cells 

Lavers 

Cells 

Base 

3 

19,278 

3 

1 , 239 

Scheme  B 

3 

17,943 

3 

1,110 

Kinetic  routines 

13.  The  water  quality  model  study  focused  on  dissolved  oxygen  resources 
and  flushing  characteristics  of  the  harbor  system.  The  WQM  simulated  the  fol¬ 
lowing  variables:  DO,  carbonaceous  biochemical  oxygen  demand  (CBOD) ,  ammonia 
nitrogen  (NH(,-N)  ,  nitrite  plus  nitrate  nitrogen  (NO2  +  NO3-N)  ,  algal  biomass 

as  carbon,  orthophosphate  (POi,-P),  and  a  conservative  tracer.  Initial  temper¬ 
atures  were  specified  horizontally  constant  yet  vertically  stratified  based  on 
the  water  quality  sampling  program  conducted  during  August  1987  by  Tekmarine , 
Inc.  However,  temperatures  varied  temporally  through  the  specification  of 
ocean  boundary  temperatures.  Salinity  was  specified  temporally  and  spatially 
constant  at  32  ppt.  The  water  quality  sampling  program  revealed  rather 
homogenous  salinities  except  near  the  Terminal  Island  Treatment  Plant  (TITP) 
effluent  and  the  Los  Angeles  River.  Salinities  were  used  in  the  WQM  for  cal¬ 
culating  DO  saturation.  Global  specification  of  a  constant  salinity  was  ade¬ 
quate  since  a  2-ppt  variation  in  salinity  at  the  temperatures  observed  results 
in  only  a  1-percent  variation  in  DO  saturation. 

14.  The  water  quality  kinetic  algorithms  were  adapted  from  the  WES  two- 
dimensional,  laterally  averaged  model  of  hydrodynamics  and  water  quality, 
CE-QUAL-W2  (Environmental  and  Hydraulics  Laboratories  1986)  and  the  HydroQual , 
Inc.,  Potomac  Eutrophication  Model  (Thomann  and  Fitzpatrick  1982).  The  water 
quality  kinetic  routines  are  detailed  in  Appendix  A. 

Transport  Calibration 

15.  Transport  calibration  consisted  of  two  studies:  comparison  of  a 
simulated  tracer  injection  using  both  the  WQM  and  CH3D,  and  a  comparison 
between  15-  and'BO-min  CH3D  time  averaging. 
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Simulated  tracer  injection 

16.  A  test  was  conducted  to  ensure  that  the  transport  properties  of 
CH3D  were  maintained  with  the  interfacing.  After  CH3D  was  calibrated,  a 
tracer  injection  was  simulated.  The  tracer  was  injected  into  the  four  surface 
cells  overlain  by  one  WQM  cell  in  the  outer  harbor  and  simulated  for  3  days. 

17.  The  same  tracer  injection  was  simulated  using  the  WQM.  For  this 
test  the  CH3D  data  were  time -averaged  over  60-min  intervals.  Plate  4  repre¬ 
sents  the  WQM  column  where  the  surface  cell  was  injected  with  tracer.  The 
circles  represent  the  arithmetic  average  of  the  six  hydrodynamic  cells  over- 
lain  by  the  single  WQM  cell.  The  figure  demonstrates  nearly  equivalent 
decline  in  tracer  in  both  the  HM  and  WQM.  The  figure  also  reveals  that  verti¬ 
cal  advection  wa  insignificant  in  mixing  the  tracer  among  layers.  The  verti¬ 
cal  diffusion  was  set  to  0.0  in  both  models  and  contributed  nothing  to 
vertical  mixing. 

Hydrodynamic  model  time  averaging 

18.  Subtidal  oscillations,  characterized  by  a  pulsating  flow  pattern 
with  areas  of  flow  direction  reversal,  were  simulated  in  CH3D.  The  flow  pul¬ 
sations  occurred  at  a  frequency  of  1  hr.  The  influence  of  the  subtidal  flow 
oscillations  on  the  HM  time -averaging  interval  was  investigated.  A  comparison 
between  15-  and  60-min  HM  time-averaging  consisted  of  two  tests.  The  first 
test  consisted  of  loading  all  WQM  cells  uniformly  with  a  tracer  at  a  concen¬ 
tration  of  10.  The  ocean  boundaries  were  0.0  concentration.  The  second  test 
consisted  of  injecting  tracer  in  the  surface,  middepth,  and  bottom  cells  of  a 
water  column  in  the  outer  harbor  and  recording  concentrations  over  24  hr.  The 
tracer  injection  corresponded  to  the  location  of  the  CH3D  and  WQM  transport 
comparison  and  was  located  in  the  "transient"  gyre  simulated  north  of  Angels 
Gate . 

19.  The  first  test  revealed  that  15-min  averaging  resulted  in 
0.2-percent  greater  mass  loss  through  the  ocean  boundaries  than  60  min  averag¬ 
ing  after  24  hr  of  simulation.  Plate  5  represents  the  WQM  cell  used  for  the 
CH3D  and  WQM  transport  comparison.  The  15-min  results  are  presented  with 
solid  lines  while  the  60-min  results  are  dotted.  The  15-  and  60-min  results 
are  nearly  equivalent.  The  conclusion  from  these  tests  was  that  60-min  hydro- 
dynamic  model  averaging  is  equivalent  to  15-min  averaging. 
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Flushing  Studies 


20.  Fxushing  studies  consisted  of  insertion  of  a  conservative  tracer 
and  noting  the  movement  and  dilution  of  the  tracer.  The  flushing  studies 
provided  a  qualitative  comparison  between  Base  and  Scheme  B  conditions.  A 
decrease  in  the  flushing  rate  prolongs  the  period  of  time  that  oxygen¬ 
demanding  substances  exert  their  influence  on  the  DO  concentration.  A 
decrease  in  flushing  rate  can  intensify  existing  water  quality  problems  and 
indicates  that  more  detailed  water  quality  analyses  are  required.  In  this 
study,  the  transport  comparisons  provided  identification  of  areas  within  the 
harbor  that  exhibited  major  decreases  in  flushing.  Areas  that  exhibited 
decreases  in  flushing  were  selected  for  more  detailed  characterization  during 
subsequent  DO  simulations. 

21.  Because  the  WQM  uses  flows  from  the  calibrated  CH3D  model,  there  is 
little  or  nothing  to  calibrate  in  the  WQM  when  simulating  a  conservative 
tracer.  The  only  parameters  that  can  be  adjusted  are  the  diffusion  coeffi¬ 
cients.  Two  field  dye  studies  were  performed  for  San  Pedro  Bay  to  verify  the 
WQM  transport;  however,  the  dye  studies  were  not  simulated  because  the  WQM 
grid  was  too  coarse  to  resolve  the  small  dye  cloud.  Field  dye  studies  did 
provide  an  order  of  magnitude  estimate  for  the  horizontal  diffusion  coeffi¬ 
cient  of  1  m^/sec  (Fischer  et  al .  1979).  Simulations  were  made  where  the  dif¬ 
fusion  coefficients  were  varied,  and  the  results  were  found  to  be  relatively 
insensitive.  This  is  not  surprising  as  many  surface-water  systems,  including 
San  Pedro  Bay,  are  advection-dominated,  with  diffusive  transport  playing  a 
minor  role.  The  vertical  diffusion  coefficient  used  was  0.1  cm^/sec .  This 
value  is  typical  of  deeper  regions  of  lakes  and  oceans  (Lerman  1971) .  The 
selection  of  this  value  represented  a  bias  toward  maximizing  DO  depression  in 
the  bottom  waters  through  inhibited  vertical  diffusive  transport. 

22.  Four  flushing  comparisons  between  Base  and  Scheme  B  are  presented: 
(a)  Tracer  Simulation  1  -  insertion  of  tracer  in  all  WQM  cells  interior  to  tlie 
breakwater,  (hi  Tracer  Simulation  2  -  insertion  of  tracer  in  the  East  Basin 
channel,  (c)  Tracer  Simulation  3  -  insertion  of  tracer  in  the  embavment  adja¬ 
cent  to  the  outer  harbor  located  between  West  Basin  of  .Midtilc'  Harbor  and  Fisii 
HarbcM'  (Seaplane  Anclu.u'age )  ,  and  (ci)  Tracer  Simulation  1  -  insertion  of  trcU’er 
in  the  West  Basin  of  Middle  liarbor  (N'avv  Harbor).  Iracer  Sinitilation  1  identi¬ 
fied  areas  ot  le..,.s  lUisiiinr,,  while  Tracer  S  i  iiiti  la  t  i  ons  throur.h  a  ex.ir,:  i  tied 

; iO  t. e  n '  ’  a  1  :  or  .i  i  w.  1 1  e  i  o'  ra  l  i  :  v  p  lo  lii  i  ,  irs  s’ ich  .i s  ic c  i  1 1.  n t  .i  1  s i >  i  i  i  s  .  i'l  r.-  ha  s  : 


Basin  Channel  comparison  was  selected  because  of  minimal  flushing  observed 
from  the  simulations;  the  Seaplane  Anchorage  was  selected  because  Scheme  B 
appears  to  isolate  these  waters;  and  the  Navy  Harbor  was  selected  because 
initial  tracer  comparisons  indicated  decreased  flushing  with  harbor 
enhancement . 

23.  The  boundaries  were  specified  exterior  to  the  breakwaters.  Tracer 
could  exit  the  outer  harbor  through  the  breakwater  openings,  but  only  water 
without  tracer  material  could  enter  the  outer  harbor  from  the  ocean  boundary. 
The  initial  tracer  concentration  was  set  to  10.0. 

24.  Plates  6-11  display  the  dilution  of  tracer  in  the  Base  condition 
while  Plates  12-17  display  the  dilution  of  tracer  in  the  Scheme  B  condition 
for  Tracer  Simulation  1,  insertion  of  tracer  in  all  WQM  cells.  Examination  of 
the  figures  reveals  that  circulation  through  Los  Angeles  Main  Channel  is 
rather  static,  but  slightly  counterclockwise  in  the  Base  condition  and  clock¬ 
wise  in  the  Scheme  B  condition.  The  clockwise  circulation  in  the  Scheme  B 
condition  is  apparent  through  the  movement  of  tracer  into  Inner  and  Middle 
Harbors.  It  should  be  noted  that  flushing  in  West  Basin  of  Middle  Harbor  is 
reduced  for  Scheme  B. 

25.  Comparison  of  the  Base  and  Scheme  B  conditions  for  Tracer  Simula¬ 
tion  2,  insertion  of  tracer  in  East  Basin  channel,  is  displayed  in 

Plates  18-23  for  the  Base  condition  and  in  Plates  24-29  for  the  Scheme  B 
condition.  Tracer  Simulation  2  corroborates  the  clockwise  circulation  in 
Scheme  B  and  the  counterclockwise  circulation  in  the  Base. 

26.  The  tracer  concentrations  in  the  five  cells  corresponding  to  the 
initial  tracer  injection  in  the  East  Basin  channel  (Plates  18  and  24)  were 
averaged  and  plotted  as  a  function  of  time  (Plate  30) .  Examination  of 
Plate  30  reveals  that  the  East  Basin  channel  flushes  more  rapidly  in  the 
Scheme  B  condition  than  in  the  existing  condition.  Tracer  concentration  in 
the  existing  condition  was  asymptotically  approaching  20  percent  of  initial 
concentration  after  25  days.  In  contrast,  tracer  surface  layer  concentrations 
achieved  20  percent  within  7  days  in  Scheme  B.  Bottom  layer  flushing  was 
slower  than  the  surface;  the  bottom  layer  required  9  days  to  achieve  20  per¬ 
cent  of  initial. 

27.  Tracer  Simulation  3,  insertion  of  tracer  between  West  Basin  of 
Middle  Harbor  and  Fish  Harbor  (Seaplane  Anchorage,  Plates  31-42),  revealed 
that,  in  the  Base,  flushing  occurred  into  the  outer  harbor  and  through  both 
Angels  and  Queens  Gates.  In  Scheme  B,  circulation  is  largely  through  Angels 
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Gate  and  into  the  Los  Angeles  Main  Channel.  The  tracer  concentrations  in  the 
six  WQM  cells  interior  to  the  embayment  were  averaged  and  plotted  as  a  func¬ 
tion  of  time  (Plate  43).  Plate  43  indicates  that  flushing  occurs  slightly 
less  rapidly  in  Scheme  B  (20  percent  in  14  days)  than  in  the  existing  condi¬ 
tion  (20  percent  in  10  days) . 

28.  The  result  of  Tracer  Simulation  4,  insertion  of  tracer  in  the  West 
Basin  of  Middle  Harbor  (Navy  Harbor),  is  presented  as  a  plot  of  concentration 
versus  time  (Plate  44) .  All  cells  in  the  Navy  Harbor  were  initialized  with  a 
concentration  of  10.0.  The  plots  were  the  average  of  the  cells.  Plate  44 
indicates  that  20  percent  initial  concentration  was  achieved  in  25  days  in 
Scheme  B  and  that  20  percent  was  achieved  in  16  days  for  the  existing  condi¬ 
tion.  The  conclusion  is  that  flushing  occurs  less  rapidly  in  Scheme  B  than  in 
the  existing  condition. 

29.  The  four  flushing  studies  indicated  areas  that  exhibited  decreased 
flushing.  Nine  locations  were  selected  for  more  detailed  characterization 
during  the  DO  simulations  (Figure  2).  The  selected  locations  are  noted  as  X-1 
through  X-9.  The  locations  prefixed  with  the  letter  I  or  B  represent  the 
interior  and  boundary  stations  sampled  during  August  1987.* 

Water  Quality  Studies 

30.  The  water  quality  studies  focused  on  dissolved  oxygen.  The  WQM 
simulated  the  following  variables;  DO,  CBOD,  NH^-N,  NO2  +  NO3-N,  algal  bio¬ 
mass  as  carbon,  and  P0i,-P.  Algal  biomass  was  reported  in  chlorophyll-a.*  The 
relationship  1  mg  carbon  =  35  mg  chlorophyll-a  (Thomann  and  Fitzpatrick 
1982)  was  used  to  relate  chlorophyll-a  to  algal  carbon.  The  kinetic  routines 
are  detailed  in  Appendix  A.  The  kinetic  constants  used  are  listed  in  Table  2. 
Boundary  conditions 

31.  Boundary  conditions  included  observed  water  quality  at  the  ocean 
boundary,*  sediment  oxygen  demand  (SOD)  at  the  bottom  (measured*),  light 
exchange  through  the  surface,  and  water  quality  of  the  Terminal  Island  Treat¬ 
ment  Plant  discharge. 


*  Tekmarine,  Inc.  1987.  "Water  Quality  Data  Collection  Program,  San  Pedro 
Bay,  California,"  report  prepared  for  US  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS,  and  US  Army  District,  Los  Angeles,  Los  Angeles,  CA. 
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Figure  2.  Water  quality  stations 


Table  2 

Kinetic  Constants 


_ Description _ 

Solar  radiation  attenuation  constant 
Fraction  solar  radiation  absorbed  at  the  surface 
Light  intensity  at  maximum  algal  photosynthesis 
Fractional  daylength 

Half -saturation  constant  for  nitrogen 

Half-saturation  constant  for  phosphorus 

Saturated  algal  growth  rate  at  20°  C 

Temperature  coefficient  for  algal  growth 

Algal  respiration  rate  at  20°  C 

Temperature  coefficient  for  algal  respiration 

Algal  maximum  excretion  rate 

Algal  mortality  rate  at  20°  C 

Algal  settling  velocity 

Phosphorus  to  carbon  ratio 

CBOD  deoxygenation  rate  at  20°  C 

Temperature  coefficient  for  CBOD  deoxygenation 

Half-saturation  constant  for  CBOD 

Ratio  of  ultimate  CBOD  to  CB0D5 

Oxygen  to  carbon  ratio 

Nitrogen  to  carbon  ratio 

Nitrification  rate  at  20°  C 

Temperature  coefficient  for  nitrification 

Half - saturation  constant  for  nitrification 

Te-  ^rature  coefficient  for  SOD 

Half-saturation  constant  for  SOD 


_ Unit _ 

0.50/m 

0.45 

20.0  W/m^ 

0.5 

0.025  g  N/m^ 

0.001  g  P/m^ 

2 . 0/day 
1.068 

0.1875/day 
1.045 
0 . 03/day 
0 . 02/day 
0 . 1  m/day 
0.025  g  PO^-P/g  C 
0 . 2/day 
1.047 

0.5  g  Oz/m^ 

1.85 

2.67  g  Oa/g  C 
0.25  g  N/g  C 
0 . 1/day 
1.08 

2.0  g  02/m^ 

1.08 

0.5  g  Oz/rn^ 


\  •/ 


32.  Solar  radiation  and  wind  shear  at  the  water  surface  were  computed 
from  weather  information.  The  computed  solar  radiation  values  were  used  in 
the  algal  growth  computations,  and  the  wind  shear  information  was  used  for 
computing  the  reaeration  coefficient  (see  Appendix  A) .  Meteorological  infor¬ 
mation  was  acquired  from  the  National  Oceanic  and  Atmospheric  Administration's 
Tape  Deck  No.  1440  WBAN  Hourly  Surface  Observations  for  the  Long  Beach  Airport 
(Station  No.  5085).  Computations  for  solar  radiation  were  based  on  the  Heat 
Exchange  Program  (Program  No.  722-F5-E1010)  available  from  the  US  Army  Engi¬ 
neer  District,  Baltimore,  Baltimore,  MD.  Code  modifications  were  necessary  to 
ensure  units  consistent  with  the  model  requirements. 

33.  The  dissolved  oxygen  at  the  ocean  boundary  was  specified  by  lin¬ 
early  interpolating  in  both  space  and  time  the  measured  data.  The  constitu¬ 
ents  algae,  P0/,-P,  NHi,-N,  NO3-N,  and  CB0D5  were  set  constant  (Table  3).  Exam¬ 
ination  of  the  water  quality  grids  (Plates  2  and  3)  indicates  that  the  ocean 
boundary  is  seaward  of  the  three  boundary  stations.  The  boundary  stations 
were  used  as  control  points  for  adjustment  of  the  boundary  values. 

Table  3 

Ocean  Boundary  Constituent  Concentrations  (g/m^) 


Variable 

Value 

Algae 

0.0175 

PO4-P 

0.0625 

NH;,-N 

0.025 

NO3-N 

0.025 

CB0D5 

0.05 

DO 

Interpolated 

34.  Examination  of  measured  SOD  data*  indicated  that  values  were  less 
than  2.0  g/m^/day.  The  exceptions  were  the  average  measured  value  of  2.1  at 
the  interior  station  in  the  East  Basin  channel  and  the  trial  estimate  of  2.1 
at  the  western  side  of  Outer  Harbor.  All  simulations  used  a  "worst-case" 
global  value  of  2.0  g/m^/day. 

35.  Constant  surface  reaeration  was  assumed.  Application  of  O'Connor's 
(1983)  relationship  between  wind  speed  and  reaeration  revealed  that  the  daily 
average  reaeration  coefficient  for  San  Pedro  Bay  varied  between 
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0.96  and  1.02  m/day;  a  constant  reaeration  coefficient  of  1.0  m/day  was  used 
for  the  water  quality  simulations. 

36.  The  Terminal  Island  Treatment  Plant  waste  effluent  was  simulated 
for  the  Base  condition.  Monitoring  reports  of  the  TITP  provided  daily  esti¬ 
mates  of  B0D5 ,  weekly  estimates  of  ammonia,  and  monthly  estimates  of  nitrate. 
Examination  of  the  ammonia  and  nitrate  values  revealed  that  the  data  were 
quite  scattered  (Table  4) .  The  nitrite  plus  nitrate  estimate  for  August  1987 
was  the  maximum  value  measured  during  the  period  September  1986  through  August 
1987. 

37.  Phosphorus  concentration  in  the  final  effluent  was  not  measured.  A 
phosphate-phosphorus  concentration  of  6  g/m3  in  the  final  effluent  from  the 
TITP  was  assumed.  This  concentration  is  typical  of  secondarily  treated  efflu¬ 
ent  (HydroQual,  Inc.  1987).  The  simulated  TITP  effluent  discharge  is  summ¬ 
arized  in  Table  5. 


Table  4 

Ammonia  and  Nitrite  +  Nitrate  TITP  Effluent  Discharge 


NH^-N,  g/sec 

NO^ 

-t-  NO^-N,  g/sec 

Date 

Averaee 

Maximum 

Minimum 

Average 

Maximum  Minimum 

September  1986- 

4.9 

18.9 

0.1 

9.3 

15.5  1.4 

August  1987 

August  1987 

0,9 

2.3 

0.2 

15.5 

Table  5 

Simulated  Effluent  PischarEe 


Variable 

CB0D5 

NH,-N 

NO2  +  NO3-N 
P0,-P 


Value .  kg/dav 
1,380 
78 
1 , 340 
430 
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Initial  conditions 


38.  Initial  conditions  were  specified  by  assuming  horizontally  constant 
yet  vertically  stratified  water  quality  based  on  the  water  quality  sampling 
program  conducted  during  August  1987  by  Tekmarine,  Inc.  The  initial  values 
represent  averages  measured  during  the  first  week  of  August  1987  (Table  6) . 

No  algae  or  CB0D5  was  detected  during  the  first  week  of  sampling.  The  monthly 
average  data  values  of  algae  were  used  as  initial  conditions.  The  initial 
phosphorus  concentration  was  substantially  inflated  by  the  values  measured 
near  the  TITP  sewage  effluent.  The  "inflated"  phosphorus  initial  concentra¬ 
tions  were  used  in  the  Scheme  B  simulations.  Temperature  was  not  simulated. 
The  temperature  values  represent  the  layer  averages  measured  during  the  month 
of  August  1987  at  the  interior  stations. 

Table  6 

Initial  Constituent  Concentrations 


Value 


Constituent 

Surface 

Middenth 

Bottom 

Temperature,  °C 

18.4 

17.2 

16.1 

Algae 

0.24 

0.08 

0.03 

P0«-P 

0.16 

0.12 

0.13 

NH^-N 

0.002 

0.002 

0.02 

NO2  +  NO3-N 

0.02 

0.03 

0.06 

CB0D5 

0.0 

0.0 

0.0 

DO 

9.4 

8.1 

6.7 

Water  quality  results 

39.  The  WQM  was  calibrated  against  observed  data  obtained  by  Tekmarine, 
Inc.,  during  August  1987.  Plates  45-90  display  simulated  algae  (Alg) ,  ortho¬ 
phosphate,  ammonia  nitrogen  nitrite  plus  nitrate  nitrogen,  biochemical  oxygen 
demand,  and  dissolved  oxygen  at  the  surface  and  bottom.  The  solid  line  repre¬ 
sents  the  Base  condition,  and  the  dashed  line  represents  Scheme  B,  The  dotted 
line  represents  the  analytical  detection  limit,  and  tiie  circles  represent  the 
observed  values.  The  station  numbering  corresponds  to  the  identification  in 
the  water  quality  sampling  program,  with  the  addition  of  nine  locatioTis  pre¬ 
fixed  with  "X"  (Figure  2).  Water  quality  station  1-3  is  not  represented  in 


Scheme  B  because  of  landfill.  Stations  I  - 10  and  I- 11  were  sampled  for  SOD  an'^ 
sediment  organic  nitrogen;  therefore,  water  column  constituents  were  not  mea¬ 
sured  and  are  not  displayed  on  the  plates. 

40.  The  measured  and  simulated  results  are  in  general  agreement  except 
for  station  1-7.  Apparently  the  Los  Angeles  River  is  contributing  some  flows 
to  the  bay  that  were  not  modeled.  The  water  quality  sampling  program  con¬ 
ducted  by  Tekmarine,  Inc.,  revealed  less  saline,  nutrient  -  enriched  surface 
water  at  station  1-7  (Figure  2). 

41.  Discrepancies  at  several  stations  between  computed  and  observed 
algae  concentrations  (such  as  1-4  in  East  Basin  channel,  Plate  48)  are  due  to 
the  use  of  "global"  algae  and  nitrogen  nutrient  values  which  exceed  initial 
measured  values  at  some  stations.  The  "global"  initial  values  represented  the 
average  over  all  interior  stations  for  each  layer.  Algae  were  nitrogen 
limited  in  the  simulations;  slight  variations  in  initial  nitrogen  concentra¬ 
tions  resulted  in  variations  in  the  simulated  algal  concentrations.  However, 
the  simulated  DO  was  insensitive  to  minor  variations  in  initial  nitrogen  con¬ 
centrations,  The  excess  nitrogen  resulted  in  rapid  growth  of  algae  followed 
by  a  gradual  decline.  Greater  algal  growth  in  the  existing  conditions,  par¬ 
ticularly  at  stations  X-6  and  X-7,  are  due  to  TITP  effluent  contributions  of 
nitrogen. 

42.  The  apparent  discrepancies  between  simulated  Base  condition  and 

measured  algal  biomass  such  as  displayed  at  the  surface  at  station  I-l 

(Plate  45)  are  not  of  concern.  Measured  algal  biomass  varies  between  0.0  and 
3  3 

0.5  g  C/m  (0  to  14  mg  Chl-a/m  ),  and  simulated  algal  biomass  was  rather  con- 

3  3 

stant  near  0.2  g  C/m  (6  mg  Chl-a/m  ).  Water  quality  standards  for  algal 

biomass  do  not  exist.  However,  it  is  generally  accepted  that  algal  concentra- 

3  3 

tions  greater  than  25  mg  Chl-a/ra  (0.875  g  C/m  )  are  undesirable.  Both  mea¬ 
sured  and  simulated  algal  concentrations  were  below  the  criterion  of  25  mg 
3 

Chl-a/m  and  much  less  than  the  general  visible  algal  concentration  of  100  mg 
Chl-a/m^ . 

43.  Examination  of  Plates  45-90  reveals  that  both  observed  and  simu¬ 
lated  DO  decreased  at  all  stations  during  August  1987.  The  decrease  in  DO  was 
due  to  a  decrease  in  boundary  DO  concentration.  For  example,  the  measured 

boundary  surface  layer  DO  .u  Angels  Gate  (station  B-1,  Figure  2)  decreased 

3  3 

from  9.0  g/m  on  4  August  1987  to  7.4  g/m  on  25  August. 

44.  Simulated  DO  of  Scheme  B  was  either  equal  to  or  less  than  existing 

3 

M.a-:iinuin  deviation.s  ot  0.5  g/m  occurred  in  the  Inner  Harbor-Bcick 

1  V 


cond i t ions . 


Channel -Middle  Harbor  of  the  Port  of  Long  Beach  (stations  1-5.  1-6,  I -11,  and 
X-5).  It  should  be  noted  that  Tracer  Simulation  1  revealed  decreased  flushing 
in  these  areas. 

45.  The  bottom  waters  exhibited  lower  DO  relative  to  the  surface 
waters.  The  maximum  deviations  between  surface  and  bottom  waters  occurred  in 
Cerritos  Channel  (station  1-10),  Back  Channel  (station  1-6),  Middle  Harbor 
(station  I-ll),  East  Basin  of  Middle  Harbor  (station  X-9) ,  and  the  dead-end 
channels  connected  to  Inner  Harbor  (stations  X-4  and  X-5).  It  should  be  noted 
that  only  the  stations  in  Back  Channel  (1-6)  and  Middle  Harbor  (I-ll) 
exhibited  differences  between  existing  and  enhancement  conditions.  Differ¬ 
ences  in  DO  between  existing  conditions  and  Scheme  B  were  not  greater  than 
3 

0.5  g/m  .  The  minimum  predicted  DO  at  all  stations  and  all  depths  was 
3  3 

6.0  g/m  .  The  minimum  of  6.0  g/m  was  predicted  in  the  bottom  layer  in  West 
Basin  of  Middle  Harbor  (station  X-7). 


Summary  and  Conclusions 


46.  The  results  of  this  model  study  indicate  that  the  Phase  1  of 
Scheme  B  for  Los  Angeles  and  Long  Beach  Ports  harbor  enhancement  will  reduce 
circulation  and  flushing  in  several  areas,  such  as  West  Basin,  Middle  Harbor, 
and  Seaplane  Anchorage.  Residual  circulation  in  the  Main  Channel  is  expected 
to  change  from  counterclockwise  (existing  conditions)  to  clockwise  for  the 
enhancement  scheme,  but  the  flushing  of  the  channel  is  noticeably  increased 
from  existing  conditions.  Scheme  B  will  result  in  less  rapid  flushing  in  the 
Inner  Harbor-Back  Channel -Middle  Harbor  of  the  Port  of  Long  Beach. 

47.  The  main  DO  impacts  of  the  enhancement  plan  would  be  experienced  in 

the  Inner  Harbor-Back  Channel -middle  Harbor  of  the  Port  of  Long  Beach.  Simu- 

3 

lated  differences  did  not  exceed  0.5  g/m  .  The  bottom  waters  exhibited  lower 
DO  than  the  surface  waters,  but  simulated  concentrations  were  greater  than 
6.0  g/m^ . 
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Appendix  A:  Kinetic  Routines 


1.  Ten  state  variables  are  included  in  the  Los  Angeles  and  Long  Beach 
Harbors  eutrophication  model: 

a.  Temperature. 

b.  Salinity. 

c.  Tracer. 

d.  Suspended  solids. 

£.  Phytoplankton. 

£.  Dissolved  inorganic  phosphorus  (PO^-P). 

£.  Ammonia  (NH^-N) . 

h.  Nitrate  +  nitrite  (NO^-N  +  NO2-N) . 

Five-day  carbonaceous  biochemical  demand. 

Dissolved  oxygen. 

2.  The  water  quality  studies  simulated  a  tracer  phytoplankton,  dis¬ 
solved  Inorganic  phosphorus,  ammonia,  nitrate  +  nitrite,  5-day  carbonaceous 
biochemical  oxygen  demand,  and  dissolved  oxygen.  Temperature  and  salinity 
were  held  constant.  Suspended  solids  was  not  modeled. 

3.  The  following  description  Includes  only  the  source/sink  terms.  The 
advectlon  and  diffusion  terms  are  omitted  In  the  state  variable  equations. 
Temperature 

4.  Temperature  kinetics  are  formulated  as  surface  heat  exchange  coupled 
with  absorption  of  solar  shortwave  radiation  (Environmental  and  Hydraulics 
Laboratories  1966*).  The  heat  exchanged  in  the  surface  layer  is  equal  to  the 
contribution  due  to  the  surface  heat  exchange  minus  the  solar  radiation  (S) 
that  passes  through  the  surface  layer.  A  fraction  (3)  of  the  solar  radiation 
that  arrives  at  the  water  surface  is  absorbed  in  the  top  approximately  0.6  m. 
The  remaining  shortwave  radiation  decays  exponentially  with  depth  according  to 
Beer’s  Law.  These  processes  are  expressed 


3TV 

3t 


A  K(E  -  T)  -  (1  -  6)  ^  exp  (-yh) 


A  Ul  -  6)  ^  exp  (-Yd)[l  -  exp  (-Yh)]| 


Surface  segment 

Otherwise 


(Al) 


*  See  References  at  the  end  of  the  main  text. 


Al 


where 


T  =  water  temperature,  °C 

V  =  segment  volume, 
t  =  time,  sec 

A  =  surface  area,  m^ 

K  =  coefficient  of  surface  heat  exchange,  m/sec 

E  =  equilibrium  temperature,  °C 

B  =  fraction  of  S  absorbed  at  the  surface 

2 

S  =  solar  shortwave  radiation,  W/m 

3 

P  =  water  density,  1,000  kg/m 
C  =  heat  capacity,  4,186  J/(kg  °C) 

Y  =  attenuation  coefficient,  1/m 
h  =  segment  height,  m 

d  =  depth  of  segment  surface,  m 
Salinity 

5.  Salinity  is  a  conservative  material  and  has  no  internal  sources  and 
sinks.  Therefore, 


(A2) 


where  f  is  salinity,  in  parts  per  thousand. 

Tracer 

6.  Tracer  is  coded  as  a  conservative  material  similar  to  salinity 
(Equation  A2) . 

Suspended  solids 

7.  Suspended  solids  is  coded  as  a  conservative  material  similar  to 
salinity  (Equation  A2) . 

Phytoplankton 

8.  One  system  is  used  to  represent  phytoplankton  assemblages.  The 
aggregation  forces  a  gross  representation  of  phytoplankton  dynamics  and 
increased  reliance  on  site-specific  characteristics  for  system  parameteriza¬ 
tion  and  interpretation.  However,  the  inclusion  of  phytoplankton  does  allow 
simulation  of  the  interactive  dynamics  of  nutrients,  phytoplankton,  and  dis¬ 
solved  oxygen.  The  rate  equation  for  phytoplankton  is  expressed 


A2 


K  =  phytoplankton  growth  rate,  1/sec 
S 

=  phytoplankton  respiration  rate,  1/sec 
=  phytoplankton  excretion  rate,  1/sec 
=  phytoplankton  mortality  rate,  1/sec 
=  phytoplankton  settling  rate,  m/sec 

9.  The  phytoplankton  growth  rate,  K  ,  is  computed  by  modifying  an 

8 

optimum  growth  rate,  G^q  ,  at  20°  C  by  the  environmental  factors  temperature, 
light,  and  available  nutrients.  The  relationship  is  given  by 

T-20 

(AL.Aj^,.Ap)  (A4) 

where 

G2Q  =  phytoplankton  maximum  growth  rate  at  20°  C,  1/sec 
0  =  temperature  coefficient  for  phytoplankton  growth 

=  light  limitation  multiplier 
=  nitrogen  limitation  multiplier 
Ap  =  phosphorus  limitation  multiplier 

10.  The  rate  multiplier  for  light  limitation  is  based  upon  the  Steele 
(1962)  function  where 

L  =  f  exp  /-  ^ 
s  \  s 

where 

L  =  limiting  factor 

2 

I  =  available  light,  W/m 

2 

1^  =  saturating  light  intensity,  W/m 

11.  The  above  expression  allows  for  simulation  of  photoinhibition  at 
light  Intensities  greater  than  the  saturating  value.  Light  penetration 
decreases  with  depth: 


+  1 


(A5) 


A3 


I  =  (1  -  3)S  exp  (-yd) 


12.  The  average  effect  of  light  on  phytoplankton  growth  In  a  particular 
segment  can  be  obtained  by  combining  Equations  A5  and  A6  and  integrating  over 
the  segment  depth  to  obtain  (Chapra  and  Reckhow  1983) 


Al  =  [exp  (-02)  -  exp  (-a^] 

yh 


where 

f  =  fractional  daylength  (0  <  f  <  1) 

(1  -  6)S  , 

j -  exp  (-yd) 

s 

“2  "  X  ~  [-Y(d  +  h)] 

s 

13.  The  rate  multiplier  for  nitrogen  and  phosphorus  limitation  is  com¬ 
puted  using  the  Monoc  relationship.  The  Monod  relationships  for  nitrogen  and 
phosphorus  are  written 


A  — 

N  ~  H  +  N 
n 


~  H  +  P 
P 


where 


=  (NO2-N  +  NO^-N)  +  NH^-N,  g/m 

=  nitrogen  half-saturation  constant,  g/m^ 

3 

=  phosphorus  concentration,  g/m 

3 

=  phosphorus  half-saturation  constant,  g/m 
.  Phytoplankton  respiration,  ,  is  written 


T-20 

K  =  R„„0 
r  20  r 


(AlO) 


where 


AA 


R2g  =  phytoplankton  respiration  rate  at  20°  C,  1/sec 
0^  =  temperature  coefficient  for  phytoplankton  respiration 

15.  Phytoplankton  excretion,  ,  is  evaluated  using  an  inverse  rela¬ 

tion  to  the  light  multiplier: 


K 

e 


(1  - 


(All) 


where  A^  is  the  maximum  excretion  rate,  1/sec.  Phytoplankton  mortality  is 
assumed  constant. 

16.  Phytoplankton  that  settle  from  a  segment  serve  as  a  source  for  the 
segment  below.  Phytoplankton  that  settle  from  a  bottom  segment  are  accumu¬ 
lated  in  a  sediment  compartment.  At  present,  however,  no  sediment  mechanism 
for  degradation  with  nutrient  release  is  planned. 

Dissolved  inorganic  phosphorus 

17.  Phosphorus  is  assumed  to  be  in  a  completely  available  form  as 
orthophosphate  (PO^)  for  uptake  by  phytoplankton.  The  sources  of  phosphorus 
are  phytoplankton  respiration  and  CBOD  mineralization;  the  phosphorus  sink  is 
phytoplankton  growth.  The  rate  equation  for  phosphorus  is  given  by 
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H  +  0 
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where 

=  ratio  of  phosphorus  to  carbon 

Jt  V 

K  =  CBOD  decay  rate  at  20°  C,  1/sec 
v20 

0^  =  temperature  coefficient  for  CBOD  decay 
0  =  dissolved  oxygen  concentration,  g/m^ 

3 

H  =0-  half-saturation  constant  for  CBOD  decay,  g/m 
V  2  3 

T  =  carbonaceous  biochemical  demand,  g  O^/m 
6c  “  ratio  of  ultimate  CBOD  to  CBODc 
6^„  =  ratio  of  oxygen  to  carbon 

VJv^ 

Ammonia 

18.  The  sources  of  ammonium  nitrogen  are  phytoplankton  respiration  and 
CBOD  mineralization.  The  sinks  of  ammonium  nitrogen  include  phytoplankton 
uptake  and  nitrification.  The  rate  relations  are  expressed 
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+  0 


where 


-  =  ammonium  nitrogen  concentration,  g/m 
=  ammonia  preference 

H  +  H 


3 


HR 


(Hn  +  B)(H^  +  n)  (E  +  n)(H^  +  n) 

K  =  nitrate  nitrogen  concentration,  g/m' 
ratio  of  nitrogen  to  carbon 
nitrification  rate  at  20°  C,  1/sec 


^20 

=  temperature  coefficient  for  nitrification 
=  0^  half-saturation  constant  for  nitrification 

Nitrate 

19.  The  system  represents  nitrite  plus  nitrate  nitrogen.  The  source  of 
nitrate  is  nitrification,  while  the  sink  is  phytoplankton  uptake.  The  rate 
equation  is  written 
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CB0D5 

20.  The  oxidation  of  carbonaceous  material  is  the  classical  BOD  reac¬ 
tion.  The  model  simulates  the  5-day  CBOD  as  the  Indicator  of  equivalent  oxy¬ 
gen  demand  for  carbonaceous  material  (HydroQual ,  Inc.  1987).  Although  the 
model  simulates  5-day  CBOD,  it  is  the  ultimate  CBOD  which  affects  the  dis¬ 
solved  oxygen  system  and  scales  nutrient  regeneration  from  the  carbonaceous 
component.  Sources  of  CB0D5  are  phytoplankton  excretion  and  mortality;  the 
sink  is  oxidation.  The  rate  equation  Is  written 
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Dissolved  oxygen 


21.  Sources  of  dissolved  oxygen  Include  surface  reaeration  and  phyto¬ 
plankton  growth.  Dissolved  oxygen  sinks  are  algal  respiration,  nitrification, 
oxidation  of  CBOD,  and  sediment  oxygen  demand. 

22;  A  distinction  is  made  between  phytoplankton  growth  using  CO^  and  NH^ 
and  phytoplankton  growth  using  CO^  and  NO^  (HydroQual,  Inc.  1987).  The  rate 
equation  is  written 
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+  A 


R(0^  -  0)  -  fij-Vo' 


where 


R  =  reaeration  coefficient,  m/sec 

3 

0^  =  dissolved  oxygen  saturation,  g/m 
0^  =  temperature  coefficient  for  SOD 

=  0^  half-saturation  constant  for  SOD  decay,  g/m^ 
Z  =  sediment  oxygen  demand  at  20°  C 


